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A selected set of topics along the borderline between perturbative and nonper- 
turbative QCD in exclusive reactions are studied. Specific problems, related 
to different mechanisms of momentum transfer to an intact hadron, are dis- 
cussed. Calculations of the space-like form factors of the pion and the nucleon 
are reviewed within a convolution scheme of short-distance (hard) and large- 
distance (soft) contributions which takes into account soft gluon emission and 
the intrinsic transverse hadron size. The failure of this scheme to reproduce 
the existing experimental data signals sizeable higher-order perturbative cor- 
rections (a K-factor of order two) and/or higher- twist contributions. 



1 Introduction 

This article presents an overwiew of exclusive processes, focusing on their end-point be- 
havior. To set the stage, we discuss and review problems relating to the (momentum) 
scales involved in form factor calculations: scale locality, infrared (IR) safety, gluonic 
radiative corrections, and the role of hadronic size effects. These issues are more pre- 
cisely described in terms of the essential mechanisms of momentum transfer to an intact 
hadron. We then use detailed calculations to investigate how these effects influence the 
predictions for F.j^[Q'^), G\j{Q'^), and G\[{Q'^) relative to existing data. 

The application of perturbative Quantum ChromoDynaniics (pQCD) to inclusive 
processes has been very successful and predictive and there is now ample experimental 
verification for a variety of reactions. In contrast, exclusive processes, though of basic 
importance for a deeper understanding of confinement, are yet not so rigorously estab- 
lished. In order that pQCD becomes applicable at the amplitude level, a short-distance 
part of the strong-interaction amplitude has to be isolated. This is then amenable 
to perturbative analysis within a hard-scattering scheme. Beyond this, however, one 
has to use additional (unsettled) nonperturbative methods to model the hadron wave 
functions which encode bound-state features. Since hadron wave functions appear in 
integrated quantities they are not directly accessible to experiment. 

Once factorization of regimes has been accomplished, renormalization group (RG) 
techniques can be employed to calculate the evolution behavior of the factorized parts. 
The logarithmic scaling violations are found to be controlled by the same nonsinglet 
anomalous dimensions as in deep- inelastic scattering 

*Invited plenary talk presented at HADRON STRUCTURE '96. High Energy Interac- 
tions: Theory and Experiments, Stara Lesna, High Tatra, Slovac Republic, February 
12-17, 1996 



1 



2 Factorization in exclusive processes 

Factorization theorems are of central importance in quantum field theory. The basic 
idea is that one can separate high-momentum from low-momentum dependence in a 
multiplicative way. For example, proving that ultraviolet (UV) divergences occuring in 
Feynman graphs can be absorbed into multiplicative renormalization factors (infinite 
constants) is instrumental in establishing renormalizability of the theory. The techni- 
cal difficulty is to prove factorization of a particular QCD process to all-orders in the 
coupling constant going beyond leading logarithms ||^. These difhculties derive from 
the fact that in QCD a new type of IR-divergence is encountered, the collinear diver- 
gence, and that in higher orders the self-coupling of gluons becomes important in the 
exponentiation of IR-divergences. 

The realization of factorization when applying to elastic form factors can be written 
in the form of a convolution of a hard-scattering amplitude (dubbed Tfj) describing the 
short-distance quark-gluon interactions, and two soft wave functions corresponding to 
the incoming and outcoming hadron [|l],Q . Generically, 



where m sets the typical virtuality in the soft parts and Q is the (external) scale charac- 
teristic of the hard (parton) subprocesses. The matching scale /i at which factorization 
has been performed is arbitrary and, assuming that /i ^ Ul, it can be safely identified 
with the renormalization scale - unavoidable in any perturbative calculation - by virtue 
of the RG equations. In this way, F can be rewritten as a function only of the coupling 
constant operative at that scale. 

As long as scaie locality is preserved, i.e., the variation of the effective coupling 
constant with U is governed by the same momentum scale, and the limit Ul — >■ 
is finite, Eq. (|l|) is valid because intrusions from the hard into the soft regime are 
prohibited. This means that Th is insensitive to long-distance interactions, i.e., it is 
IR safe. All IR-sensitivity resides in the hadron distribution amplitudes ^^"(o*'*) which 
are independent of large momentum scales and may depend on the external scale 
only through RG-evolution. Both the subtraction procedure of UV poles in the soft 
parts and the cancellation of IR divergences in the hard part are not uniquely fixed. 
Nevertheless, they have to ensure that the asymptotic behavior of F is IR-insensitive 
and governed by the leading anomalous dimensions associated with vertex and quark 
self-energy corrections. 

Adopting a factorization scheme, the initial (final) state of the hadron has a certain 
probability distribution for finding its valence quarks carrying longitudinal momentum 
fractions < = / < 1 in a — *• OO frame. Apart from the slow per- 
turbative Q'^-evolution, this randomness depends only on the uncalculable confinement 
dynamics and not on the specific hard-scattering collision, i.e., it is universal. Hence 
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where = 2, C = 1 for the meson and iV = 3, C = 2/3 for the nucleon, respectively, 
and 7_p is the anomafous dimension associated with quark self-energy. 

To solve the evolution equation, for hadron H has to be expressed as an 

orthogonal expansion in terms of appropriate hypergeometric functions which constitute 
an eigenfunction basis of the gluon-cxchange potential, i.e., 

Q') = $if (X.) £ 5f Ha:.) exp I /f %A9m)\ , 

^0 IV /i J 

(3) 

where is the asymptotic amplitude (see below) over fractional momenta propor- 

tional to the weight w{Xi) of the orthogonal basis, and ^^^^ denote the corresponding 

eigenfunctions. The coefhcients -B^^-* of this expansion are associated with matrix el- 
ements of composite lowest-twist operators with definite anomalous dimensions taken 
between the vacuum and the external hadron. In the meson case (leading twist=2), the 
eigenfunctions of the diagonalized evolution equation are the Gegenbauer polynomi- 
als C^^^ 1^, which correspond to conformal operators |^ with associated anomalous 
dimensions given by 
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Introducing the relative coordinate ^ = Xi — orthogonality with respect to the 
weig ht w{i) = (1 - i^) yields expansion coefficients going like 



In 



A2 



2(2n-f 3) 
(2 + n)(l + n) 7-1 



diC'f{i)^^-\i,Q'), (5) 



meaning that for higher orders they decrease like Xjv? , provided ^^'^•'(xi, /i) < K.x\ 
as Xj ^ for some e > ||l|. 

In the nucleon case (leading twist=3), an orthogonal normalized basis of the evo- 
lution kernel is provided by linear combinations of Appell polynomials which depend 
on two variables j^J^ . Here orthogonality alone is insufficient to fix the eigenfunctions 
uniquely. It was first shown in Q that the expansion coefhcients Bn become analyt- 
ically tractable up to any desired polynomial order M > i -|- J in terms of strict 
moments 



[dx] x\ X3 $Ar {Xk , jJ, 



of the mixed-symmetry nucleon distribution amplitude ||] $ 
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The projection coefficients Cl^-, the anomalous dimensions of trihncar twist-three quark 
operators 

7f ^ = \ (^^^ + 2^"^^) 

(where are the zeros of the characteristic polynomial that diagonalizes the evolution 
equation) , and the normalization constants A^^ up to order M = 4 are tabulated in . 

The moment values - calculated, for instance, via QCD sum rules [^,0 - provide 
only local constraints that are not sufficient for the distribution amplitude to be re- 
constructed in an unambigous way | [lf| ] . One can always add some oscillating function 
which vanishes at the points fixed by the local constraints but which contributes out- 
side. Thus one has to impose global constraints as well, for shaping the distribution 
amplitude as a wJiote. Such constraints have been successfully used in Refs. - ||l^ to 
ensure dominance of the lowest-order contributions and minimize the influence of dis- 
regarded higher-order terms. In a complete set of nucleon distribution amplitudes 
was determined which satisfy existing QCD sum rules pO| , p^ with comparable degree 
of accuracy while avoiding unphysical oscillations. These solutions organize themselves 
across a "fiducial orbit" that is characterized by a scaling relation between the form- 
factor ratio I I / and the expansion coefficient ^4 (cf. Eq. (|)). 

Hadron distribution amplitudes derived this way from QCD sum rules show an 
asymmetric balance of longitudinal momentum fractions of valence quarks. Convoluted 
with the corresponding hard-scattering amplitudes they lead to form factors which have 
the right magnitude and QCD-evolution behavior Q . On the other hand, it was pointed 
out in 18 1 that asymmetric distribution amplitudes enhance the contributions of end- 
point regions and that, extracting these regions, the leading perturbative contribution 
to the form factor is reduced to a small fraction. This depletion of the form factor indi- 
cates sensitivity to the gluon offshellness in the end-point region, presumed to be large. 
Hence the perturbative treatment turns out to be inconsistent, meaning that uncalcu- 
lated higher-order terms may be important, even leading. To reinstate the validity of 
pQCD, end-point contributions have clearly to be suppressed. To this end, a modified 
convolution scheme [|9| ~ still within the hard-scattering picture - will be discussed 
below which incorporates Sudakov suppression due to radiative gluon corrections. 



3 Momentum transfer mechanisms 

There are basically two schemes for describing the transfer of a large external momentum 
Q to an intact hadron during elastic scattering: hard-gluon exchange and the Feynman 
mechanism. We begin with the first one which is tightly connected to pQCD and relies 
on the factorization theorem. Following this rationale, the struck quark connects to the 
other valence quarks via highly off-shell gluon propagators, meaning that the transverse 
interquark distances are rather small, viz. of order \/ Q and that all partons share 
comparable fractions of longitudinal momentum. Thus Th can be reliably computed 
as a power series in the running coupling constant [0J|] (see the Ihs of Fig. |l|) . 

In the asymptotic limit Q'^ — > OO, only the n = term in Eq. (||) survives so that 
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the RG-asymptote of the distribution amphtude reads (70 < 7n for aU ?2 > 0) 

<^>if(Q'-oo) = 5rMx.) lim In-^ (9) 

Q^^oo y 1\QCD J 

with C = Cp (pion) and C = 3Cp/2 — 2Cb (nucleon), where Bq^^ is the hadron 
wave function at the origin of coordinate space and the hmit of the logarithm amounts 
to the wave- function renormahzation factor Z2- From these distribution amphtudes 
one infers that asymptotically the most likely configurations are those in which the 
valence quarks share longitudinal momentum in a uniform way, i.e., Xi = 1/2 for the 
pion and Xj = 1/3 for the nucleon. Within this scheme, when confinement sets in, 
a quark is not able to venture too far from the antiquark (in the pion) or the other 
two valence quarks (in the nucleon). This poses constraints on the offshellness of the 

involved propagators, typified by XiX'j > (^^qcd/Q'^^ ^xp (47r//?as). The problem 
then is that in the end-point regions Xi,x'j — > 0, 1 the gluons become nearly on-shell 
(i.e., real) and the hard-gluon exchange mechanism becomes unreliable. This is also 
reflected in the behavior of form factors. As a rule, narrow distribution amplitudes 
yield for reasonable values of Ctg results which are unrealistically low to be consistent 
with the data. Obviously, broad distribution amplitudes are required to match the data. 

One method |^ extracts distribution amplitude moments from QCD sum rules 
using local vacuum condensates. In the pion case, the large values obtained this way 
for the moments can only be realized by a "double-humped" distribution amplitude of 

the form 4>cz{^) — '^^fnXil — a;)(l — 2x) which is end-point dominated, whereas 
the central region (Xj = 1/2) is depleted. This distribution amplitude yields a pion 
form factor in much better agreement with data but at the expense that it accumulates 
its main contributions from the end-point region where a perturbative treatment is less 
reliable. 

The approach taken in [ pO| is conceptually quite different. To avoid the inherent 
deficiencies of moment inversion (see for criticism p8| ) , the pion distribution amplitude 
was computed directly from QCD sum rules via dispersion relations. This approach 
makes use of nonlocal vacuum condensates which afford for the finite average virtuality 
of the vacuum quarks. Parametrizing the nonlocal quark condensate by a Gaussian, a 
model distribution amplitude was obtained which gives lowest moments < > with 
n = 2, 4, 6 close to those of ^^jj = /^6x(l — x), but which has a significantly wider 
shape and no dip in the central region. 

Concerning the nucleon the situation is technically more complicated, though, per- 
haps, less prone to criticism. Up to now, all derived model distribution amplitudes (see, 
e.g., @,|l^) rely on moment-inversion techniques within rather large uncertainty inter- 
vals. However, a consistent pattern has emerg ed which seems to encapsulate the 
main characteristics of the true nucleon distribution amplitude. 

Let us consider now the other basic mechanism for elastic scattering due to Feyn- 
man pl|. In this scheme, almost all of the hadron's momentum is carried off by a 
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Fi gure 1: Mechanisms for momentum transfer during elastic scattering. The Ihs 
shows hard-gluon exchange within pQCD. The blob S containing soft gluon lines (and 
an analogous one with soft quark-antiquark lines not shown here) spoils factorization 
but is power-suppressed, i.e., non- leading. The (rhs) shows the Feynman mechanism 
using for purposes of illustration quark and gluon lines. The leading quark is denoted 
by a heavy line, while all other lines represent wee quarks and soft gluons. 

single parton, the others being "wee" . This picture is consistent with a configuration in 
which only the struck quark is within an impact distance 1/Q of the electron while all 
other partons have rather random positions in the transverse direction, building a soft 
"cloud" with transverse size ^ \/Q. Once the elastic scattering has happened, rear- 
rangements are necessary to change quarks and gluons into hadrons. This conversion 
procedure (visualized on the rhs of Fig. |l|) is controlled by the overlap of the initial and 
final state wave functions and cannot be computed within pQCD. 

For the remainder of this report we will consider only calculations which are based 
on the hard-scattering picture. 

4 Modified Convolution Scheme 

In deriving Eq. (|l|) , we tacitly assumed that the A;_|_-dependence of the quark and gluon 
propagators in Tjy can be ignored. This is tantamount to factorizing the /c_|_-dependence 
into the distribution amplitudes which are the wave functions integrated over up 
to the factorization scale. Then, in the limit — OO, the only gluon radiative cor- 
rections remaining uncancelled are those giving rise to wave-function renormalization. 
Hoewever, in the end-point region the parton transverse momenta in Tjj cannot be a 

priori ignored since, say, for the pion, -\- ^ X^x'jQ'^. As a result, the 

transverse distance between the quark and the antiquark becomes large compared to 
1 /Q and the corresponding gluon line is no more part of the hard-scattering process but 
should be counted to its soft part. In other words, the hard-gluon exchange mechanism 
should be replaced by that of Feynman. 

The physical basis of the modified convolution scheme (MCS) jl^ is to dissect the 
process in such a way, so that for transverse distances large compared to 1/Q (play- 
ground of the hard-scattering mechanism) but still small relative to the true confinement 
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regime - characterized by 1 / ^QC D ~ the hadron wave function is modified to exhibit 
the effect of Sudakov enhancements expHcitly up to the transverse scale retained in Th. 
Going over to transverse configuration space, the modified wave function reads 

^'^,A^^^/KQ,^iren) = e-^^^''\x,,l/h) , (10) 



where the factor 



N 



exp {-S) = 



i=l 



s{xi,bi,Q) + / _ —iqigiJi )) 



+ Xi^ x\ (111 



comprises gluon corrections in terms of the functions s(Xj, 6j, Q) |T9j and accounts 
for RG-evolution from the IR-scale l/fej to the renormahzation scale [iren via the 
quark anomalous dimension (in the axial gauge) ^q[g(jj'^)) = —ds/ll -\- 0(«^). 
The explicit expressions for the Sudakov functions are given in [p2| . The Sudakov 
exponential factor resums contributions from two-particle reducible diagrams (giving 
rise to double logarithms), whereas two-particle irreducible diagrams (giving rise to 
single logarithms) are absorbed into the hard scattering amplitude Tjj ||2^. It can be 
conceived of as a Unite renormahzation factor to the hadron wave function The 
leading double logarithms derive from those momentum regions where soft gluons (all 
four-momentum components small) and coUinear gluons to the external quark lines 
overlap. These contributions arc numerically dominated by the term 



2Cf 6Q , In {i,QlV2KQCD) 
exp < — In^ In — —. — r^, > (12) 



where is one of the fractions Xj or X^, and (3 = (33 — 2nj)/3 is the first-order term of 
the Gell-Mann and Low function encountered before. The single logarithm stems from 
the running coupling constant and the double logarithm contains the exponentiated 
higher-order corrections - required by RG - rendered finite by the inherent IR-cutofff 
l/6j. This marks the crucial difference between the MCS and previous approaches 
dealing with isolated quarks where such IR-cutoff parameters had to be introduced as 
external regulators. 

For small transverse distances (or equivalently, l/6j ^ ^iQ)^ gluonic radiative cor- 
rections are treated as being part of Th and are excluded from the Sudakov form factor. 
Consequently, for < \/2/biQ the Sudakov functions s{^i,bi, Q) are set equal to 

zero. On the other hand, as 6j increases decreases, reaching zero at biAgcD = 1- 
In the pion case, there is only one transverse scale, notably, the quark- antiquark sepa- 
ration b, and suppression is automatically accomplished. Indeed, when it happens that 
one Sudakov function s(^C,,bi, Q) = (or equivalently that the corresponding expo- 
nential is set equal to unity) the other (negative) Sudakov function in the exponent, 
s(l — ^, b, Q), diverges, thus providing sufficient suppression. 
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Our approach [gij to the choice of the appropriate IR-cutoff in calculating observ- 
ables with several transverse momentum scales involved - like nucleon form factors - 
is to postulate that long wave-length (compared to the typical interquark separations) 
gluons "see" the nucleon as a whole, i.e., in a color-singlet state and cannot resolve its 
color structure. This is technically implemented by setting b = max{6i, 62; ^3} = 

61 = 62 = ^3 ("MAX" prescription). As a result, (i) logarithmic a^-singularities 
in the end-point region are screened by the exponentially decreasing Sudakov factors, 
(ii) the form-factor integrands are IR-safe for all possible kinematic configurations and 
receive perturbative contributions which saturate, i.e., which are rather insensitive to 
distances of order 1/Aqcd- It was outlined in (2^Jl^ that other choices may lead to 
uncompensated «s-singularities. In what follows we present those aspects of the cal- 
culation which differ from the standard case. We also present theoretical predictions 
in comparison with available experimental data. A more detailed level of description is 
given in |^,|l5|. 

All told, the nucleon form factor recast in the transverse configuration space reads 



1 fi 

Gm(Q^) = — / [dx][dx' 



3 Jo 



X exp 



(47r)2 {Any 



-Sj{^i,bi,Q,fi) 



J2 Tjix, x', b, Q, /i)y,(x, x', b, ji) 
j 

(13) 



where the transverse separation vectors between quarks 1 and 3, 2 and 3, and 1 and 2 

are defined as follows: 61 = b'^, 62 = ^2' ^3 ~ ^2 ~ ^1- The modified wave functions 
are given by {Xi = 1 — Xi) 

^Pin) ^ _l_|4(_2)^-3^123 + 4(-2)^-2^132 + ^231^231 + ^r21^321 



XiX 

± 2(11 



^231^132 + ^^32^231 + ^321^123 + ^m^32l] } , 



(14) 



p(n) 



1(2) 

T'lX' 

^{2(1) 

.1"-, I- 



{3(0)^^32^132 T ^231^231 T ^231^132 T ^t32^23l} 

^321^321 + ^321^123 + ^f23^32l]±^f23^123} , (15) 



X3X1 



where the lower signs and the numbers in parentheses refer to the neutron. The diagrams 
of hard-gluon exchanges in the MCS can be conveniently combined fia] to give 
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xix'y/^' 



QbAKJ{x2x'^Y'^Qb2 



(16) 
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(17) 



where Kq is the modified Bessel function of order (the Macdonald function). The argu- 
ments of the running coupUng constant, t^j, are defined as the maximum scale of either 
the longitudinal momentum (X Q or the inverse transverse separation OC 1/bi, appear- 
ing in the argument of Kq. They are associated with the virtualities of the exchanged 

gluons, namely, t 
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max 



{xix[y/^Q, 1/6 



21 



max 



[xix[y/^Q,i/b 



and ti2 = ^22 = max 
proton wave function reads 

^123 (x, 6, /i) 



[X2x',y/^Q, 1/62 



. Note that the Fourier transform of the 



fN{fJ')^123{x, fl)ni23{x,b) 



(18) 



where the /i;_L-dependent part is modeled by 

1 



^i23ix,b) = (Att) exp 



4a2 



(xiXsbl + X2X^b\ + XiX2^\ 



(19) 



Similar, albeit simplified, expressions are obtained also for the pion form factor 
The theoretical predictions are shown in Figs. ^ and |[ 
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Figure 2: Spacelike pion form factor comprising Sudakov corrections and the intrinsic 
transverse size of the pion wave function p^ ] in comparison with experimental data. The 
curves correspond to the Chernyak-Zhitnisky model 1^ (solid line) and the asymptotic 
wave function (dotted line). The data are from |^^. 
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Figure 3: Spacelike magnetic form factors of the nucleon: proton (Ihs) and neu- 
tron (rhs) {Aqcd = 180 MeV). The open circles indicate Ff data The calcu- 
lated curves have been obtained in |2j] with the "MAX" prescription, and including 
QCD-evolution and the intrinsic fc_L-dependence of the nucleon wave function (the lat- 
ter normalized to unity). The shaded area contains the predictions derived for the set 
of nucleon distribution amplitudes determined in via QCD sum rules. 



5 Conclusions 



The above discussion raises a troubling question. If the self-consistent calculations 
within the MCS are insufficient to describe the existing data, is all hope lost for the 
applicability of the pQCD paradigm to exclusive reactions? At one level, no: higher- 
order corrections, i.e., a K-factor of order 2 - in principle, also computable within the 
hard-scattering scheme - might bridge the gap to the data. But on a deeper level, 
the news is sobering. The results may be interpreted as evidence for the failure of the 
hard-scattering mechanism in exclusive reactions at accessible momentum transfer and 
this calls for nonperturbative mechanisms. 
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